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Abstract — This paper presents several new construction tech- 
niques for low-density parity-check (LDPC) and systematic 
repeat-accumulate (RA) codes. Based on specific classes of 
combinatorial designs, the improved code design focuses on 
high-rate structured codes with constant column weights 3 and 
higher. The proposed codes are efficiently encodable and exhibit 
good structural properties. Experimental results on decoding 
performance with the sum-product algorithm show that the 
novel codes offer substantial practical application potential, for 
instance, in high-speed applications in magnetic recording and 
optical communications channels. 

Index Terms — Low-density parity-check (LDPC) code, system- 
atic repeat-accumulate (RA) code, accumulator design, sum- 
product algorithm, combinatorial design. 



I. Introduction 

STRUCTURED low-density parity-check (LDPC) codes 
have been designed recently based on finite geome- 
tries and, more generally, balanced incomplete block designs 
(BlBDs), see, for example, |lT|-||5l. One of the main ad- 
vantages of these structured LDPC codes is that they can 
lend themselves to very low-complexity encoding, as opposed 
to random-like LDPC codes. Moreover, experimental results 
show that the proposed codes perform well with iterative 
decoding. 

In the first part of this paper, we present novel high- 
rate structured LDPC codes with constant column weights 3 
and higher, based on cyclic BIBDs (CBIBDs), resolvable 
BIBDs (RBIBDs), and cycUcally resolvable cyclic BIBDs 
(CRCBIBDs). We obtain several infinite classes of (fc,r)- 
regular LDPC codes with values of k varying from 3 to 8 
(Theorems [T] m andlU as well as one infinite class of (fc,r)- 
regular LDPC codes for k any prime power (Theorem |3), all 
admitting flexible choices of r and the code length. The pre- 
sented results are more general than previous ones in [E]-||ll. 
Our proposed LDPC codes have good structural properties, 
their Tanner graphs are free of short cycles of length 4 and 
have girth 6. Their code rate is high, achieving 9/10 or higher 
already at small block lengths. This is typically the rate of 
interest in high-speed applications, for example, in magnetic 
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recording and optical communications channels. Most classes 
of the constructed codes have a quasi-cyclic (or similar) struc- 
ture, which allows linear-complexity encoding with simple 
feedback shift registers. In addition, many classes exhibit very 
sparse parity-check matrices. Experimental results on decoding 
performance show that the novel LDPC codes perform very 
well with the sum-product algorithm (J6l, significantly better 
than known BIBD-LDPC codes and random Gallager LDPC 
codes |7| for short to moderate block lengths (cf. Fig. |2}. 
We observe furthermore, very interestingly, for the proposed 
codes a performance gain as the column weights grow larger, 
and a particular good performance of LDPC codes based on 
CRCBIBDs. 

In the second part of the paper, we apply our combinatorial 
construction techniques to systematic repeat-accumulate (RA) 
codes. RA codes were first introduced in ||8] as a serial 
concatenation of a rate- ^ repetition code and a rate-1 convolu- 
tional code with transfer function yT-jy , called the accumulator 
(Fig. [TJ. Between these two constituent codes, an interleaver 
permutes the output of the repetition code and in some cases, 
the encoding scheme is also refined by a rate-a combiner 
This combiner performs modulo-2 addition on sets of a bits. 
Here, we consider systematic RA codes (abbreviated from 
now on as sRA codes), where the message of length K is 
concatenated with the output of the accumulator of length M. 
For a detailed description of the encoding process, we refer the 
reader to |9 1. The scheme leads to a parity-check matrix of the 
form H — [H1H2], where Hi is an MxK matrix with column 
weight q and row weight a, and H2 is an M x M double 
diagonal matrix designated solely by the accumulator [9|. It 
is obvious that H is sparse and thus the code can be decoded 
with the sum-product algorithm in the same way as for LDPC 
codes. 

One drawback of sRA codes is that they possess many 
weight-2 columns in H2, which has a negative effect on 
the minimal distance of the code and thus on the decoding 
performance. These low-weight columns result solely from the 
accumulator structure and hence Johnson and Weller [10] have 
introduced a new accumulator type to increase the weights of 
these columns. The new accumulator produces many columns 
of weight 3 instead of weight 2 and, consequently, the resulting 
codes are called weight-3 RA (w3RA) codes. In the present 
paper, we propose a generalized accumulator structure, leading 
to novel coding schemes, termed weight-q RA (wqRA) codes. 
This allows code constructions for higher column weights. 
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Fig. 1: Low-complexity encoder of systematic RA codes 



having the advantage that their parity-check matrices are even 
more similar to those of regular LDPC codes and thus have 
a comparable gain in decoding performance (cf. Fig. |5). We 
note that, while BIBDs by their axioms can always be used to 
build up parity-check matrices of sRA codes, it is the specific 
nature of CBIBDs, RBIBDs and CRCBIBDs that is exploited 
to build up parity-check matrices of wqRA codes. 

The paper is organized as follows: Section |II] briefly in- 
troduces BIBDs and their usefulness in general LDPC code 
design. We give furthermore a short overview of previous 
results on LDPC code design based on CBIBDs, RBIBDs 
and CRCBIBDs. In Section [nil we extend these previous 
results and present several novel high-rate LDPC codes with 
higher column weights. Experimental results on decoding 
performance are given in Subsection IIII-DI In Section |IV] we 
propose a generalized accumulator structure for sRA codes, 
leading to our novel coding schemes of wqRA codes. We 
develop several code constructions from CBIBDs, RBIBDs 
and CRCBIBDs that are capable of producing such wqRA 
codes. Simulation results are presented in Subsection IIV-DI 
The paper is concluded in Section |V] 

II. General LDPC Code Design Based on Balanced 
Incomplete Block Designs 

A. Balanced Incomplete Block Designs 

We give some standard material on combinatorial designs 
that is important for our further purposes. For encyclopedic 
references, we refer the reader to (ITTI . (IT2l . Various applica- 
tions of combinatorial designs to coding, cryptography and 
information security are described in llT3l , l]T4|, including 
many further references. 

Let X be a set of v elements and B a collection of 
fc-subsets of X. The elements of X and B are called points 
and blocks, respectively. An ordered pair {X,B) is defined 
to be a balanced incomplete block design (or 2-design), and 
denoted by BIBD(ti, k, A), if each pair of points is contained 
in exactly A blocks. A BIBD(w, A;, 1) is also called a Steiner 
2-design (or system). It is straightforward that in a BIBD each 
point is contained in the same number r of blocks, and for the 
total number h of blocks, the parameters of a BIBD satisfy 
the relations bk = vr and X{v — 1) = r{k — 1). For {X, B) a 
BIBD(w, k, A), its incidence matrix is avxb matrix A = {a,ij) 
[l < i < V, I < j < b), in which a^ = 1 if the i-th point of 
X is contained in the j-th block of B, and a^ = otherwise. 
Clearly, A is unique up to column and row permutation. 

Let {X, B) be a BIBD(ti, A;, A), and let u be a permutation 
on X. For a block B = {bi,...,bk} G B, define B'' := 
{bl, . . . , bl}. If B" := {B'' : B C,B} = B, then a is called 
an automorphism of {X,B). If there exists an automorphism 



a of order v, then the BIBD is called cyclic and is denoted by 
CBIBD(w, fc, A). In this case, the point-set X can be identified 
with I^y, the set of integers modulo v, and a can be represented 
by cr : i —> i + 1 (mod v). For a block B — {5i, . . . , 6^} in a 
CBIBD(w, k, A), the set B + i := {bi + i (mod v),...,bk + i 
(mod v)} for i e Z„ is called a translate of B, and the set 
of all distinct translates of B is called the orbit containing 
B. If the length of an orbit is v, then the orbit is said to be 
full, otherwise short. A block chosen arbitrarily from an orbit 
is called a base block (or starter block). If k divides v, then 
the orbit containing the block B = {O, f , 2f , . . . , (fc - l)f } 
is called a regular short orbit. For a CBIBD(i;, k,l) to exist, 
a necessary condition is w = 1 or A (mod k{k — 1)). When 
1) = 1 (mod k(k — 1)) all orbits are full, whereas if v = k 
(mod k{k — 1)) one orbit is the regular short orbit and the 
remaining orbits are full. 

A BIBD is said to be resolvable, and denoted by 
RBIBD(t;, k, A), if the block-set B can be partitioned into 
classes TZi , ■ . ■ , TZr such that every point of X is contained 
in exactly one block of each class. The classes TZi are called 
resolution (or parallel) classes. If TZi is a resolution class, 
define TZf := {B'^ : B e 7?.,}. An RBIBD is called 
cyclically resolvable if it has a non-trivial automorphism a 
of order v that preserves its resolution {TZi, . . .TZr}, i.e., 
{TZl,...TZ'^} = {TZi,...TZr} holds. If, in addition, the 
design is cycUc with respect to the same automorphism a, 
then it is called cyclically resolvable cyclic, and denoted by 
CRCBIBD(w, k, A). 

In a CBIBD(w, k, 1), we can define a multiset AB := {bi — 
bj : i,j = l,...,k;i ^ j} of differences for a base block 
B = {bl, . . . , bk}. Let {Bi}i^j, for some index set /, be all the 
base blocks of full orbits. If w = 1 (mod fc(fc— 1)), then clearly 
IJjgj ABi = Z,„ \ {0}. The family of base blocks {Bi}i^j 
is then called a (cyclic) difference family in Zy, denoted by 
CDF(w,fc,l). 

Let k be an odd positive integer, and p= I (mod fc(fc — 1)) 
a prime. A CDF{p, k, 1) is said to be radical, and denoted by 
RDF(p, k, 1), if each base block is a coset of the fc-th roots 
of unity in Zp [TSl . 

B. BIBD LDPC Code Design 

A (binary) LDPC code |7| is a linear block code with a 
sparse parity-check matrix (over GF{2)). If the parity-check 
matrix has both constant column and constant row weights, 
then the code is said to be regular, otherwise irregular. We 
summarize the state-of-the-art of general LDPC code design 
from BIBDs with respect to our further purposes (cf. LI |-[5l, 
|16|): The v x b incidence matrix of a BlBD{v,k, 1) defines 
a parity-check matrix H of a binary {k, r — 6A:/t;)-regular 



TABLE I: Possible exceptions: An RBIBD(t;, fc, 1) with k 
V = k (mod k(k — 1)). 
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BIBD-LDPC code of length b = {v{v - l))/{k{k - 1)), if we 
identify the v points with the parity-check equations and the h 
blocks with the code bits. In this case, all column weights of H 
are equal to k, and all row weights are equal to r. The dimen- 
sion of the code is K ~ 6— rank2(i?), where rank2(i?) denotes 
the 2-rank of H, and the code rate is i? = (6 — rank2(i/))/6. 
In general, often the calculation of rank2 [H] depends on the 
specific structure of the BIBD(w,fc,l). It is clear from the 
classical Fisher's Inequality that rank2(il) < v. Furthermore, 
it is elementary to see that, if 2 | |5t' ^^^ rank2(if ) > w— 1, 
with equahty if and only if 2 | fc (cf. fTT, Thm. 2.4.1]). More 
precise results are in particular known for BIBD (i;, 3, l)s 1 18l. 
As in a BIBD(u, fc, 1) no pair of points have more than 
one block in common (i.e., no pair of columns of H contains 
more than one "1" at the same positions), the Tanner graph of 
a BIBD-LDPC code is free of cycles of length 4 and hence 
has girth at least 6. Moreover, the fact that every pair of points 
occurs in exactly one block ensures that the girth is exactly 
6. Since each column in H has weight fc and no two columns 
have more than one "1" in common, there are fc parity-check 
equations that are orthogonal on every code bit (i.e., the code 
is one-step majority logic decodable). Thus, the minimum 
distance of the code is at least fc + 1. The same lower bound 
holds for the size of minimal stopping sets I.19J . 

Some specific classes of BIBDs proved to be especially 
useful in the code design: 

(a) LDPC code constructions based on CDFs have been con- 
sidered in |2)-|4| for constant column weights fc = 3 to 
7. CDF-LDPC codes have the structural advantage that 
each orbit of a base block in the BIBD corresponds to a 
circulant submatrix in the parity-check matrix. Hence, 
these codes a have quasi-cyclic structure (cf. Il20l ). 
which allows linear-complexity encoding with simple 
feedback shift registers. 

(b) LDPC codes based on RBIBDs have been constructed 
in H, |J5l for constant column weights fc = 3,4. 
RBIBD-LDPC codes have the property that they have 
very sparse parity-check matrices. 

(c) LDPC codes based on CRCBIBDs have been designed 
in |4|, |5| for constant column weights fc = 3,4. 
CRCBIBD-LDPC codes exhibit very sparse parity-check 
matrices and, in addition, their structure allows linear- 



complexity encoding with simple feedback shift registers 
similarly as for quasi-cyclic codes. 

III. New LDPC Codes with Higher Column Weights 

In this section, we extend the previous results on LDPC 
codes based on CDFs, RBIBDs and CRCBIBDs to higher 
column weights. Experimental results on decoding perfor- 
mance with the sum-product algorithm show that the novel 
LDPC codes offer substantial practical application potential 
(Subsection UlLDl) . 

A. CDF LDPC Codes 

Theorem 1: Let p be a prime. Then there exists a {k,r)- 
regular quasi-cyclic LDPC code of length pr/k and rate at 
least (r — fc)/r based on a CDF(p, fc, 1) for the following 
cases: 

(1) {k,p) = (3, 6i + 1) for any positive integer t, 

(2) {k,p) ~ (4, 12i + 1) for any positive integer t, 

(3) (fc,p) = (5, 20i + 1) for any positive integer t, 

(4) {k,p) = (6, 30i + 1) for any positive integer t 7^ 2, 

(5) {k,p) — (7, 42i+ 1) for any positive integer t > I with 
the possible exceptions p = 127,211 as well as primes 
p e [261239791, 1.236597 x lO^^] such that (-3)^4^ = 
1 in Zp. 

(fc,p) — {8,p) for all values of p = 1 (mod 56) < 10^ 
with the possible exceptions p = 113, 169, 281, 337, 
{k,p) = (9,p) for all values of p ee 1 (mod 72) < 10^ 
with the possible exceptions p — 289, 361. 

Proof: The proof relies on known infinite ((l)-(5)) and 
finite ((6)-(7)) families of CDF(p,fc, l)'s (cf. (El] and the 
references therein; |12]). ■ 

We note that Cases (l)-(5) have been presented in JSJ-llll. 

Example 1: A (5, 26)-regular quasi-cyclic LDPC code of 
length 546 and rate at least 0.81 can be obtained from a 
CDF(105,5,1). 

B. RBIBD LDPC Codes 

Theorem 2: Let w be a positive integer. Then there exists 
a (fc, r)-regular LDPC code of length vr/k and rate at least 
{r—k)/r based on an RBIBD(u, fc, 1) for the following cases: 
(1) (fc, v) — (3, 6t + 3) for any positive integer t. 



(6) 



(7) 



TABLE II: Existence of a CRCBIBD(pfc, k, 1) with k = 5, p < 10^, and A: = 7 or 9, p < 10"* 
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(2) (fc, v) ~ (4, 12t + 4) for any positive integer t, 

(3) (fc, v) — (5, 20i + 5) for any positive integer t with the 
possible exceptions given in Table H] 

(4) (fc, v) — (8, 56i + 8) for any positive integer t with the 
possible exceptions given in Table |T] 

Proof: The proof is based on known infinite series of 
RBIBD(v, fc, l)'s (cf. [|12J . Ii22j and the references therein). 



We remark that Cases (1) and (2) of have been considered 
in|4|, |5|. 

Example 2: A (5, 46) -regular LDPC code of length 1702 
and rate at least 0.89 can be constructed from an 
RBIBD(185,5,1). 

Theorem 3: If v and k are both powers of the same prime, 
then a (fc,r)-regular LDPC code of length vr/k and rate at 
least {r — k)/r based on an RBIBD(w, k, 1) exists if and only 
if (i; - 1) = (mod {k - 1)) and w = (mod k). 

Proof: It has been shown in ll23l that, for v and k both 
powers of the same prime, the necessary conditions for the 
existence of an RBIBD(t7, k, A) are sufficient. ■ 

C. CRCBIBD LDPC Codes 

Theorem 4: Let p be a prime. Then there exists a {k,r)- 
regular LDPC code of length pr and rate at least (r — k)/r 
based on a CRCBIBD(pfc, k, 1) for the following cases: 

(1) {k,p) ~ (3, 6i + 1) for any positive integer t, 

(2) {k,p) = (4, 12i + 1) for any odd positive integer t, 
(3a) {k,p) = (5, 20i+ 1) for any positive integer t such that 

p < 10"^, and furthermore 
(3b) {k,p) = (5, 20t + l) for any positive integer t satisfying 
the following condition: e + 1 is not a 2'^+^-th power in 
Zp, or equivalently (ll + 5\/5)/2 is not a 2'^+^-th power 
in Zp, where 2^ is the largest power of 2 dividing t and 
e a 5-th primitive root of unity in Zp, 
(4) (fc,p) = (7, 42i + 1) for any positive integer t satisfying 
the following condition: there exists an integer / such 



that 3^ divides t and e + l,e'^ + e+l, 



are 



3^-th 



powers but not 3^+^-th powers in Zp, where e is a 7-th 
primitive root of unity in Zp, 
(5) {k,p) = (9,p) for the values of p = 1 (mod 72) < 10"* 

given in Table HIl 
Moreover, there exists a (fc, r)-regular LDPC code of length 
£r and rate at least (r - k)/r based on a CRCBIBD(£fc, /c, 1) 
for the following cases: 



(6) {k, t) for fc = 3, 5, 7, or 9, and £ is a product of primes 
of the form p = 1 (mod k{k—l)) as in the cases above, 

(7) {k, tj = (4, t) and ^ is a product of primes of the form 
p = 12t+l with t odd. 

Proof: The proof is given in Appendix |A] ■ 

We note that Cases (1) and (2) have been treated in f4\, 

[:5|. 

Example 3: A (5, 51)-regular LDPC code of length 2091 
and rate at least 0.90 can be obtained from a CRCBIBD (41 x 

5,5,1). 

D. Simulation Results 

For experimental results on decoding performance, we em- 
ployed the iterative sum-product algorithm, as proposed in |6|, 
on an AWGN channel with a maximum of 50 iterations per 
codeword. Fig. |2] shows the bit-error rate (BER) performance 
of LDPC codes that have been constructed using the differ- 
ent combinatorial techniques presented in this section, and 
random regular Gallager LDPC codes ||7j. A legend displays 
the following information in the respective order: code type, 
construction method in brackets, and a quadruple [N, K, R, k] 
consisting of block length N, dimension K, code rate i? and 
column weight k of the parity-check matrix. 

The left plot of Fig. |2] shows, that codes with higher column 
weights may have a relatively small performance deficit at 
low signal-to-noise-ratios (SNRs), but outperform their coun- 
terparts for higher SNRs significantly. The CBIBD-LDPC code 
based on Netto's construction (cf. i2)-|i4l) and the RBIBD- 
LDPC code of Johnson and Weller (cf. PI, fSl) have been 
truncated in length to yield parameters comparable to those 
of our new RBIBD-LDPC code. Both plots demonstrate that 
the new structured codes outperform known random Gallager 
LDPC codes and known structured BIBD-LDPC codes with 
comparable parameter sets. We observe a performance gain 
as the column weights grow larger and a particular good 
performance of CRCBIBD-LDPC codes. 

IV. New WEiGHT-q RA Codes 

As an extension of the accumulator type for w3RA codes 
introduced by Johnson and Weller [101, we present a general- 
ized accumulator for wqRA codes. This accumulator structure 
was similarly introduced by Liva et al. |24| with focus on the 
design of irregular RA codes. The intention of their work is to 
allow a more flexible choice of the nodes degree distribution 
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Fig. 2: Decoding performance of various novel structured LDPC codes 



for optimizing the code in the sense of density evolution f25l, 
while we use the accumulator to provide new sRA codes with 
an improved decoding performance closer to those of regular 
LDPC codes. Let gi, . . . , gq-i £ {1, •••, M} denote the design 
parameters of the generalized accumulator, where gi ^ gj for 
i ^ j, and s; :== Y.]=i 5j. ^ ^ {1, . . . , q - 1}. Let r^ be the 
output of the combiner at time i. The output of the accumulator 
Pi is now 



Pi 
Pi 



Ti ® Pi 



l,...,Sl 

.Si + 1,...,S2 



Pi '■ — T^i (SPi — si ffl ■•■ (S Pi 



Sq-1 



1,...,M. 



This scheme again results in a parity-check matrix of the 
form H — [H1H2], where the i-th column of H2 has 1 -entries 
at rows i and i + si, 1 < I < q~l, given that these rows exist. 
Thus, the design parameters gj specify the vertical distance 
between the j-th and {j + l)-th 1 -entry in the columns of 
H2- The i?i-part is unaffected by the accumulator and hence 
remains as introduced in Section |I] Consequently, most of the 
columns of H2 have the same (high) column weight q as the 
columns of Hi, leading to an improved decoding performance 
compared to sRA codes. 

A. CBIBD sRA Codes 

The incidence matrix of a CBIBD has a quasi-cyclic struc- 
ture A = [^1^2, ..., A;] of circulant submatrices Ai and 
can easily be transformed into the parity-check matrix of 
an sRA code. For this purpose, we consider the underlying 



difference family. A CDF F always contains a base block 
Bt G F with 1 e ABt- The translates of Bt give rise 
to the circulant submatrix At, where every column has two 
consecutive 1 -entries. We retain all pairs of ones and delete 
the others to obtain the required double diagonal matrix H2 
for sRA codes. Finally, we can use an arbitrary number of the 
remaining Ai, 1 < i < I, i ^ t to form Hi (Fig. [3a] i. 

It is also possible to construct weight-g RA codes from 
CBIBDs, where one of the design parameters gi, ..., gq-i can 
be chosen arbitrarily, say gi w.l.o.g. We have to identify the 
base block Bg with gi G AB^ in the underlying CDF and 
transform the associated As into the accumulator matrix H2 
by simply deleting the 1 -entries of the upper triangular. The 
other design parameters arise from the remaining differences 
of B,. 

1} Netto's Construction for column weights 3: Let GF{p) 
be a Galois field of prime order p = 6t + 1 and G its 
multiplicative group. Let w be a primitive element of the 
field. Then the set F := {uj' G^* mod p \ 1 < i < t} is a 
CDF(6i+ 1,3,1) Il26l. 

2) New sRA and w3RA Codes: With Netto's construction 
we can obtain both sRA codes and w3RA codes. For this pur- 



pose, we have to find the base block Bg = {0 



.s+2t |^s+4t 



*} 



(mod p), s e {1, ...,t}, that generates the differences 1 and gi, 
respectively. The position of Bg depends on the distribution 
of the primitive element w and it can be shown that this 
search problem is equivalent to solving the discrete logarithm 
(see Appendix |B]). Thus, finding the required base block in 
general is difficult for large difference families. We may now 
proceed to construct the desired sRA codes and w3RA codes 
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(a) CBIBD(19, 3, 1) sRA code (b) RBIBD(21, 3, 1) sRA code 

Fig. 3: Parity-check matrices of two systematic RA codes 



as described eariier in this section. 

3) Buratti's Construction for column weights 4 and 5: 
Let GFijj) be a Galois field of prime order p — 12t + 1 
and H its multiplicative group. If we find some b e GF{p), 
B = {0, 1, b, &2} such that A+B := {&, - bj\l < j < i < 4} 
is a representative system of the cosets of H^ in H, then 
F := {lo'^'B mod p | 1 < i < i} is a CDF(p,4,l) ||27l . 
Similarly, for p — 20t + 1 we can construct a CDF(p, 5,1) by 
F := {uj^^'B mod p | 1 < i < t}, if we find some b e GF{p), 
B = {0, 1, b, b^, b^} such that A+i? is a representative system 



nm- 



of the cosets of H^^ in 

4) New sRA, w4RA and w5RA Codes: For obtaining sRA 
codes or wqRA codes with gi = 1, we have to identify the 
base block Bi E F with 1 £ ABi. For Buratti's construction 
we can show that the desired block is the last one of the family 
given hy i — t. 

Proof: It holds that uj^p^^^/^ mod p = p — 1. In case 
of p — 12t + 1, the last block contains and w^* mod p = 
p — 1. These two elements generate the required difference 
— (p — 1) mod p = 1. In case of p = 20i + 1, the proof is 
analogous. ■ 

Therefore, we can efficiently derive sRA codes or wqRA 
codes with gi — 1 from the underlying difference family in 
contrast to the construction in Subsection IIV-A2I 

B. RBIBD sRA Codes 

Let B{i) be the column of a matrix which contains consec- 
utive ones at position i and i + 1. For sRA codes, we have to 
identify all B{i), 1 < i < A/ — 1, in the incidence matrix 
of an RBIBD for constructing the double diagonal matrix 
i?2- Unfortunately, these columns are spread over different 
resolution classes and thus, taking these columns for H2 would 
destroy several classes. Therefore, we propose to build up H2 
only from a few resolution classes and let the others unaffected 
(Fig. [3b] i. The advantage of this approach is that we can use 
an arbitrary number of the intact resolution classes to form 
Hi and hence have more flexibility in the code design. It is 
important to omit only complete resolution classes, because we 
then remove a definite number of ones per row which ensures 
the regularity of the resulting parity-check matrix. 

1) sRA Codes by Johnson and Welter: Johnson and Weller 
rely on the construction of Ray-Chaudhuri and Wilson [28 1 for 
generating an RBIBD(t;, 3, 1), also called a Kirkman triple 
system (KTS). The last three resolution classes of the KTS 
have a special structure consisting of circulant submatrices. 
These classes can be transformed into a double diagonal matrix 



for sRA codes with simple row and column permutations [TOl. 
More precisely, let 

jd 

r 
If 



c 



ja 
I" 

F 
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be the last three resolution classes, where every class is 
vertically divided into three circulant submatrices /*, which 
are cyclic shifts by i places of the A//3 x M/3 identity 
matrix. Johnson and Weller first replace I'^,!'^ and /'* with 
all zeros matrices and then permute it into a double diagonal 
matrix, where the row permutations must be applied to the 
entire incidence matrix. Finally, an arbitrary number of the 
remaining resolution classes can be used to form Hi. 

2) Modification for New w3RA Codes: We propose here not 
to replace the submatrices J"^, /^^ and I^ with zero matrices, 
but to keep the entries in C. With the same permutations we 
obtain the parity-check matrix of a w3RA code, where gi = 1 
and (72 arise from the construction. 

C. CRCBIBD sRA Codes 

Codes from CRCBIBDs combine the resolvability of 
RBIBDs and the cyclic character of CBIBDs, so that we 
gain the specific advantages of both code classes. The cyclic 
structure enables linear-complexity encoding, while we have 
the flexibility of RBIBDs to adjust the length and rate of 
the code independently by removing an arbitrary number of 
resolution classes. Although the resulting CRCBIBD codes are 
not quasi-cyclic, we may realize linear-complexity encoding 
either by a modified scheme of feedback shift registers as 
proposed in [4| or by employing the encoder of an sRA code. 
For the latter, there is no need to maintain the quasi-cyclic 
structure of the parity-check matrix and we thus may use an 
arbitrary selection of the resolution classes for designing the 
code. 

1) Construction of Genma, Mishima and Jimbo: With the 
construction of Genma, Mishima and Jimbo [29|, we can 
obtain a CRCBIBD(pfc, A:, 1) with p prime and k odd, if there 
exists an RDF(p, fc, 1). In their proof, they construct all base 
blocks of the design and use them to partition the design into 
resolution classes. For our purpose to construct sRA codes, 
we need one of the resolution classes given by 7?.' q and its 
class orbit 

^i,/ = T^'j.o + ^ (mod pk), I elk, j e I, 
where / is a suitable subset of {0, 1, ..., 2mn — 1}. 
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Fig. 4: Transformation of C into iJj for a CRCBIBD(39, 3, 1) code 



2) New sRA and wqRA Codes: The class orbit TZ', i,l G l^k 
consists of k resolution classes, which are cyclic shifts of each 
other. Let C be the incidence matrix of this class orbit. Now 
let {xi,yi) be the position of an arbitrary 1 -entry in C and 
S be the vertical distance to the next 1 -entry in column j/i 
(Fig. m. It must hold that gcd{6,pk) = 1, otherwise we have 
to choose another 1 -entry with different S. Furthermore let 11 
be a permutation vector for the rows of C, where 11 [i] is the 
value at position i, 1 < i < pk. Now we compute 

n[z] = xi + S{i — 1) (mod pk) 

and apply the permutation to the rows of the entire incidence 
matrix. Here we consider only the partial matrix H2 corre- 
sponding to C (Fig. m. For H2 it holds that 

1) there exists exactly one column B{i) with 1 -entries at 
row i and i + I (mod pk) for every i, 1 < i < pk, and 

2) all B{i) are distinct. 

Proof: The column yi of C with 1 -entries at rows xi and 
xi + 5 gives the column -8(1) of H2- Because of the cyclic 
resolvability there must be a resolution class with a column 
yt+i that contains 1-entries at rows xi + tS (mod pk) and 
xi+S + t6 (mod pk), t e 1, ...,pk - 1. With gcd{6,pk) = 1, 
we have S'Zpk — Zpj. and hence we never use a row twice. 
The columns yt+i result in B{t + 1) of H2 and thus we have 
B{i)' s for every i. Moreover, B{i) must be unique, because 
there can only exist one column with 1-entries at rows xi +r5 
(mod pk) and xi + {r + \)5 (mod pk), r e {0, ...,pk — 1}, 
due to the axioms of the design. 

The columns y, must be distinct because the blocks of TV, q 
can never have two pairs of entries with difference 5 due 
to constructional reasons. Thus, the resulting B{i) must be 
distinct. ■ 

Finally, we can transform H2 into a double diagonal matrix 
by simple column reordering and deleting the superfluous 
ones. We can even obtain wqRA codes with an arbitrary design 
parameter gi with the restriction that gcd{gi, pk) — 1. For this, 
we modify our computation to 



nf^ 



ffij 



xi + 5{i — 1) (mod pfc), 1 < i < pk. 



D. Simulation Results 

In this section, we compare the decoding performance of 
sRA codes, wqRA codes and regular LDPC codes for our 
new constructions (Fig. |5]l. These codes differ solely in the 
i/2-part of their parity-check matrices and thus, the largest 
differences in performance are expected for short to moderate 
block lengths or for lower code rates, as in these cases the 
impact of H2 is relatively higher than that of Hi. For the 
decoding we use the sum-product algorithm |6| on an AWGN 
channel with a maximum of 50 iterations per codeword. The 
LDPC codes are again described by a quadruple [N, K, R, k], 
and the sRA/wqRA codes by a slightly different quadruple 
[N,K,R,q], where q is the column weight of Hi. The first 
plot of Fig. [^demonstrates the performance gain of a new low- 
rate w3RA code at relatively small block length, compared 
to the sRA and LDPC code of the same parameters. In the 
second plot, our w5RA-code shows a similarly well low-SNR 
performance as the corresponding sRA code and the excellent 
high-SNR performance as the LDPC code. The sRA code 
suffers a high error-floor at approximately 3.75 dB, arising 
from many weight-2 columns in the parity-check matrix. This 
error-floor turns out stronger compared to the first plot, since 
we reduced the column weights of the _ff2-part from 5 (instead 
of 3) to 2 in order to obtain a double diagonal matrix. 
Consequently, the larger decrease of the column weights 
negatively affects the decoding performance in the error-floor 
region. The third plot demonstrates the performance gain of 
our novel w3RA codes compared to the sRA codes of Johnson 
and Weller, for various codes of rate 0.85 and block lengths 
varying from 1220 to 3020. The fourth plot shows a particular 
good high-SNR performance of our w5RA-code compared to 
the corresponding sRA and LDPC code. 

V. Conclusion 

We have designed in this paper new classes of structured 
LDPC codes with high code rate and low-complexity encod- 
ing. Based on specific classes of BIBDs, we obtained several 
infinite classes of (fc, r)-regular LDPC codes with values of k 
varying from 3 to 8 as well as an infinite class of {k,r)- 
regular LDPC codes for k any prime power, all admitting 
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-e — CBIBD-LDPC (Netto) (842, 421 , 0.5, 3] 

-k— CBIBD-sRA (Subsection IV-A) [842, 421, 0.5, 3] 

-0 — CBIBD-W3RA (Subsection IV-A) [842, 421 , 0.5, 3] 




-e— CBIBD-LDPC (Buratti) [2005, 1604, 0.8, 5] 



-4f— CBIBD-sRA (Subsection IV-A) [2005, 1604, 0.8, 5] 
-$ — CBIBD-W5RA (Subsection IV-A) [2005, 1604, 0.8, 5] 



2 3 4 5 6 7 

Signal-to-noise ratio E^/N^j (dB) 




2.5 3 3.5 

Signal-to-noise ratio E^/N^j (dB) 



- e - RBIBD-sRA (Johnson, Weller) [1220, 1037, 0.85, 3] 
— B — RBIBD-W3RA (Subsection IV-B) [1220, 1037, 0.85, 3] 

.|qO - m - RBIBD-sRA (Johnson, Weller) [1940, 1649, 0.85, 3] 
RBIBD-W3RA (Subsection IV-B) [1940, 1649, 0.85, 3] 

- • - RBIBD-RA (Johnson, Weller) [3020, 2567, 0.85, 3] 
— $ — RBIBD-W3RA (Subsection IV-B) [3020, 2567, 0.85, 3] 
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-e — CRCBIBD-LDPC (Genma) [1 025, 820, 0.8, 5] 
-k— CRCBIBD-sRA (Subsection IV-C) [1025, 820, 0.8, 5] 
-$ — CRCBIBD-W5RA (Subsection IV-C) [1025, 820, 0.8, 5] 
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Fig. 5: Decoding performance of various novel systematic RA and wqRA codes 



flexible choices of r and the code length. We have furthermore 
addressed sRA codes, proposing a generalized accumulator 
structure for higher column weights that replaces the con- 
ventional accumulator and leads to an encoding scheme, we 
have termed weight-g RA code. This allowed an improved 
decoding performance closer to those of regular LDPC codes, 
along with the low encoding complexity of turbo-like codes. 
Compared to sRA codes, the relatively high error floors could 



be lowered significantly. The encoding scheme is applicable to 
our new construction techniques and may be used for further 
combinatorial constructions that lack an efficient encoding. 
The presented constructions arise from cyclic and resolvable 
BIBDs and allows to use an arbitrary number of block orbits 
or resolution classes for more flexibility in the code design. 
Therefore, we can adjust the rate and length of the codes 
independently and thus produce a wider range of codes. The 



proposed novel LDPC and systematic RA codes in this paper 
offer good structural properties, perform very well with the 
sum-product algorithm, and are suitable for applications, e.g., 
in high-speed applications in magnetic recording and optical 
communications channels. 

Appendix A 
Proof of Theorem|4] 

Proof: Let p be a prime. We first assume that k is odd. 
Then, by a result of Genma, Mishima and Jimbo ||29l , a 
CRCBIBD(pfc, fc, 1) exists, whenever there is an RDF(p, fc, 1). 
The following infinite ((i)-(iii)) and finite ((iv)) families of rad- 
ical difference families are known (cf. ifTSl and the references 
therein; fT2l): 
(i) An RDF(p, 3, 1) exists for all primes p=l (mod 6). 
(ii) Let p = 20i + 1 be a prime, let 2^ be the largest power 
of 2 dividing t and let e be a 5-th primitive root of unity 
in Zp. Then an RDF(p, 5, 1) exists if and only if e + 1 is 
not a 2'^+^-th power in Zp, or equivalently (ll + 5\/5)/2 
is not a 2'^+^-th power in Zp. 
(iii) Let p = 42i + 1 be a prime and let e be a 7-th primitive 
root of unity in Zp. Then an RDF(p, 7, 1) exists if and 
only if there exists an integer / such that 3-'^ divides t 
and e + l,e^ + e + l, ^ ^_^i'^ are 3'^-th powers but not 
3''^"''^-th powers in Zp. 
(iv) An RDF(p, 9, 1) exists for all primes p < 10"' displayed 
in Table M 
Moreover, in [29] a recursive construction is given that 
implies the existence of a CRCBIBD(^fc, k, 1) whenever £ is 
a product of primes of the form p = 1 (mod k{k — 1)). In 
addition, a CRCBIBD(5p, 5, 1) has been shown i30l to exist 
for any prime p=l (mod 20) < 10"^. 

We now consider the case when k is even: In fTP\, a 
CRCBIBD(4p, 4, 1) is constructed for any prime p = 20t + 1, 
where t is an odd positive integer Furthermore, via the above 
recursive construction, a CRCBIBD(4£, 4, 1) exists whenever 
£ is a product of primes of the form p = 12t + 1 and t is odd. 
The result now follows. ■ 



Appendix B 
Search Problem 



Let Bi = {u 



,i+2t 



;'+■**} be the i-\h block of the dif- 



ference family F := {oj^G^* mod p \ 1 < « < i} over GF{p) 
with prime order p = 6i + 1. From the construction arises the 
following matrix 



Ai^ = 






UJ ' 
..c+t+1 



to 

UJ 



c+t-1 
2t-l 



W 



c+5t 



UJ 



c+5t+l 



UJ 



h6t-l 



where the i-th column consists of the differences ABi, and 
uj'^ = w(ci;^* — 1) (mod p). Note, that such an element c 
must exist, since w is a primitive element of GF{p). Let 
s G {0, ..., i— 1} denote the index of the column that contains 
entry 1 . The index s indicates the position of the desired block 
Bs containing the difference 1. Now, the following problems 
are equivalent: 



1) Computation of the discrete logarithm c = Log^{uj'^) 

2) Computation of s, such that 1 S ABg- 

Proof: First, assume that the discrete logarithm is solved 
and thus c is known. Then, s = t — c (mod t). Conversely, 
assume that s is known. Let r e {0, ..., 5} be the row index of 
the 1-entry in AF, such that AF(r,s) = 1- Since there are only 
six possible rows, we can find r in constant time. The solution 
of the discrete logarithm is then given hy c~Qt~ rt ~ s. ■ 
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